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ABSTRACT 



Context. Successful inference from asteroseismology relies on at least two factors: that the oscillations in the stars have amplitudes 
large enough to be clearly observable, and that the oscillations themselves be stable enough to enable precise measurements of mode 
frequencies and other parameters. Solar-like p modes are damped by convection, and hence the stability of the modes depends on the 
lifetime. 

Aims. We seek a simple scaling relation between the mean lifetime of the most prominent solar-like p modes in stars, and the 
fundamental stellar parameters. 

Methods. We base our search for a relation on the use of stellar equilibrium and pulsation computations of a grid of stellar models, 
and the first asteroseismic results on lifetimes of main-sequence, sub-giant and red-giant stars. 

Results. We find that the mean lifetimes of all three classes of solar-like stars scale like T^^ (where T^g is the effective temperature). 
When this relation is combined with the well-known scaling relation of Kjeldsen & Bedding for mode amplitudes observed in narrow- 
band intensity observations, we obtain the unexpected result that the height (the maximum power spectral density) of mode peaks in 
the frequency power spectrum scales as g"^ (where g is the surface gravity). As it is the mode height (and not the amplitude) that 
fixes the S/N at which the modes can be measured, and as g changes only slowly along the main sequence, this suggests that stars 
cooler than the Sun might be as good targets for asteroseismology as their hotter counterparts. When observations are instead made 
in Doppler velocity, our results imply that mode height does increase with increasing effective temperature. 

Key words, convection - stars: oscillations - data analysis 



1. Introduction 

Asteroseismology is now opening new windows on the interi- 
ors of solar-like stars, and allowing inferences to be made on 
the masses, sizes, and ages to levels of precision that would not 
otherwise be possible. These results have come from ground- 
based campaigns (Bedding & Kjeldsen 2006) and from space- 
craft observations, e.g., WIRE (Bruntt 2007), MOST (Matthews 
et al. 2007) and SMEI (Tarrant et al. 2007). CoRoT (Baglin et al. 
2006) is now producing exciting results on hot solar-like main- 
sequence stars (Michel et al. 2008); while NASA's Kepler mis- 
sion (Basri et al. 2005) will increase by more than two orders 
of magnitude the number of solar-like stars we can observe by 
seismology (Christensen-Dalsgaard et al. 2007). 

A first detailed attempt to predict amplitudes and hence the 
detectabiUty of solar-Uke p-mode oscillations in stars other than 
the Sun was presented by Christensen-Dalsgaard & Frandsen 
(1983). They made computations for a grid of stellar models 
in the colour-magnitude diagram, and concluded that both ve- 
locity and intensity oscillation amplitudes generally increased 
with increasing age and with increasing mass along the main se- 
quence. Kjeldsen & Bedding (1995) used the same model data as 
Christensen-Dalsgaard & Frandsen, and found that the p-mode 
velocity amplitudes in the models scaled as Vqsc °^ (L/MY where 
i = 1. Kjeldsen & Bedding (1995) also used linear theory and 



observational data to obtain a relation between the velocity and 
intensity amplitudes, i.e., (dL/L),i oc Vosc/T^g (assuming the am- 
plitudes in intensity were infeiTed from narrow-band observa- 
tions). 

The scaling relations of Kjeldsen & Bedding (1995) have 
been used extensively within the field of asteroseismology. 
However, several studies have since suggested modifications to 
the relations. Houdek et al. (1999) calculated a new grid of stel- 
lar models, and included the perturbations to the convective heat 
and momentum fluxes in the mode stability computations. From 
this improved grid of stellar models Houdek et al. (1999) still ob- 
tained the same scaling relation as Kjeldsen & Bedding (1995) 
for L/M less than three times the solar value, but found that 
5 = 1.29 for M/Mq > 1.4. Samadi et al. (2005), on the other 
hand, used numerical 3D simulations of near-surface convection 
to obtain s w 0.7, a result which has been supported to some 
extent by a number of ground-based observations of oscillations 
in solar-like stars. Very recently, Michel et al. (2008) concluded 
from the first solar-like oscillation observations made by CoRoT 
that intensity amplitudes were about 25% lower than the predic- 
tions of Samadi et al. (2005). 

In this letter, we continue the work initiated by Christensen- 
Dalsgaard & Frandsen (1983) and Kjeldsen & Bedding (1995). 
However, rather than looking at the mode amplitudes we try in- 
stead to find a simple scaling relation to describe the lifetimes. 
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T, of the solar-like oscillations. The solar-like oscillations are 
damped (and excited) by convection, and so the lifetimes are 
potentially an extremely useful diagnostic of near-surface con- 
vection in stars. Our first motivation in seeking a robust scaling 
relation is therefore to help elucidate the underlying physics at 
play. Our second motivation is to enable simple, but robust, life- 
time predictions to be made for main-sequence, sub-giant and 
red-giant stars. As explained below, the lifetimes affect the de- 
tectability of the modes, so that lifetime predictions have a role 
to play in target selection. 

For a given mode amplitude, strong damping spreads the 
power of the mode over a large range in frequency, implying a re- 
duction in the height of its peak in a power frequency spectrum 
compared to a mode with weaker damping. With t defined as 
the e-folding lifetime for the amplitude, the mode peak observed 
in the frequency power spectrum will have a FWHM linewidth 
of A = 1/{7tt). If a peak is well resolved, its height H may 
be expressed in terms of the mode amplitude, A (either Vqsc or 
{6L/L)^) and mode linewidth. A, according to 



Hence, when the mode is resolved, it is H that determines the 
S/N in power, not alone. 

We note that a robust predictor is also desirable to fix the 
underlying parameters of artificial asteroseismology data. Tests 
with artificial data (e.g., the hare-and-hounds exercises con- 
ducted in support of CoRoT and Kepler) are an important part 
of developing and validating data analysis codes. 

The layout of the rest of our paper is as follows. We begin in 
Section|2]with a description of the theoretically computed mode 
lifetimes, and measured mode lifetimes, that were used to infer 
the scaling relation. The measured mode lifetimes come from 
observations of 12 main-sequence, sub-giant and red-giant stars, 
made by ground-based telescopes and the CoRoT, WIRE and 
SMEI spacecraft. We then proceed in Section[3]to seek a relation 
between mode lifetimes and the fundamental stellar parameters. 
Our main result is that if we average hfetimes over the most 
prominent modes, the average lifetime scales to a good approx- 
imation as (t) oc T^^. We finish in Section |4] with a discussion 
of our scaling result, and its implications for future selection of 
asteroseismic targets. 

2. Data 

We consider mode lifetimes computed from the grid of stellar 
models used by Chaplin et al. (2008) with measured mode life- 
times from 12 solar-like stars, sub-giants and giants. 

The grid of stellar models used by Chaplin et al. had masses 
in the range 0.7 to 1 .3 Mq, and ages in the range from the ZAMS 
to 9 Gyr Padova isochrones (Bonatto et al. 2004; Girardi et al. 
2002, 2004) were used to specify the primary characteristics of 
each model, i.e., mass M, radius R, effective temperature T^^fi, 
and luminosity L. The composition was fixed at X = 0.7 and 
Z - 0.019 for all models. In order to estimate mode lifetimes for 
each model Chaplin et al. (2008) then performed stellar equi- 
librium and pulsation computations, as described in Balmforth 
(1992), Houdek et al. (1999), and ChapHn et al. (2005). The 
pulsation computations required estimates of M, Teff , L and the 
composition as input. The computations gave as output estimates 
of the powers and lifetimes of the radial p modes of each stellar 
model. We used the powers and lifetimes to compute the implied 
heights of the modes, and then estimated an average lifetime 



(t) for each model by averaging lifetimes (and linewidths) over 
the five most prominent radial modes, i.e., those with the largest 
heights. More details may be found in Chaplin et al. (2008). 

We also use measurements of mode lifetimes of solar-like os- 
cillations from observations of 12 main-sequence, sub-giant and 
red-giant stars. The measured lifetimes are shown in Table [T] 
together with estimates of stellar masses, luminosities, and ef- 
fective temperatures, the measured large frequency separations 
of the p-mode spectra, and the list of references from which the 
data were taken. 

It is important to point out that the mode hfetimes were not 
measured in the same way on all 12 stars, and the quoted av- 
erages are not entirely self consistent; nor do the quoted val- 
ues necessarily correspond precisely to an average over the five 
strongest radial modes (the datum chosen for the model data). 
However, given the level of precision in the measurements, and 
given that, in many cases, it was only possible to estimate the 
lifetimes of a small number of modes in the vicinity of the high- 
est amplitude, we feel this is not a major cause for concern. 

The mismatch between the measured lifetimes of the dif- 
ferent stars will no doubt introduce some extra scatter into the 
lifetime-stellar parameter relations we seek to constrain. We also 
note that for the stars observed by the Kjeldsen & Bedding group 
(v Indi, a Cen B, yS Hydri, Procyon and t Ceti), lifetimes were 
evaluated from frequency scatter of mode frequency ridges in the 
Echelle diagram without including the effects of rotation, which 
are known to artificially lower the value of the mode lifetime 
(see Bazot et al. 2007 for discussion). We comment in Section|3] 
below on the impact of this bias on our results. 

3. Results 

Following Kjeldsen & Bedding (1995), we first sought to obtain 
a power-law relation between the average lifetimes and L/M. 
The left-hand panel of Fig.[T]plots, as a function of L/M, the av- 
erage lifetimes from the stellar equilibrium and pulsation com- 
putations (diamond symbols), and the average lifetimes from 
the observations (crosses with error bars). The average lifetimes 
from the pulsation computations - which we recall do not in- 
clude evolved stars - follow a well-constrained power-law rela- 
tion of the form 

<t) cc (L/M)-0-5i±''''5 . (1) 

The solid line shows the best-fitting power law. The average life- 
times of the observed main-sequence stars show a rough match 
to this power law; however, the more evolved stars with large 
L/M depart significantly from the prediction. We conclude that 
a power law in L/M fails to provide an adequate scaling relation 
for main-sequence, sub-giant and giant stars. 

It struck us that in order to get the observed lifetimes to fol- 
low a single curve, we needed an independent variable which 
reversed the order on the abscissa of the main-sequence points 
and giant points. Use of T^ff as the independent variable would 
fulfill this requirement. The right-hand panel of Fig. [T]therefore 
plots the lifetime data as a function of Tgff . The pulsation com- 
putations again follow a well-constrained power law, described 
by 

(T> cx T;^-'^'-\ (2) 

but this time the observational data on the main-sequence stars, 
sub-giants and giants also show a good match to the best-fitting 
curve. We add that the model prediction for ^ Hydrae presented 
by Houdek & Gough (2002) is not inconsistent with this power 
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Table 1. Mode lifetimes and fundamental stellar parameters for 12 solar-like stars, sub-giants and giants. 



Name 


M/Mq 


L/Lq 


7'eff 


Ay 


T 


^ Hydrae 


3.07" 


61.1 + 6.2° 


5000 ± 100 K" 


7.11+6.14 yuHz" 


2 days" 


V Indi 


0.847 ± 0.043'' 


6.21 ± 0.23'' 


5291 ± 34A:'' 


24.25 ± 0.25 yuHz'' 


16.2_^ J days" 


a Cen A 


1.105 ± 0.007'' 


1.522 ± 0.03'' 


5810 ± 50 K'' 


106.2-' 


3.9+ 1.4 days« 


a Cen B 


0.934 ± 0.006'' 


0.503 ± 0.02'' 


5260 ± 50 K" 


161.38 ± 0.06 yuHz* 


3.3+';^ days'" 


j3 Hydri 


1.07 ± 0.03' 


3.51 ± 0.09' 


5872 ± 44 K' 


57.24 ± 0.16 yuHz-' 


2.3*5 j days' 


HD49933 


1.2* 


0.53+0.01* 


6780 + 130K* 


85.9 ± 0.15 yuHz* 


2.3+^-] days* 


Procyon 


1.48' 


6.9' 


6500 K' 


55.90 ± 0.08 yuHz' 


l.Stli days'" 


T Ceti 


0.783 ± 0.012" 


0.488 ± 0.01" 


5264 K" 


169 ± 0.2yuHz" 


1.7 + 0.5 days" 


eOph 


2.4 ± 0.4'' 


59'' 


4887 + 100 K" 


5.3 + 0.1 yuHz« 


2.7+551 'lays" 


Arcturus 


0.8 ± 0.3'' 


200 + 10' 


4290 ± 30 K"^ 


0.8 ^lllz' 


24.1 days'" 


y8 Ursae Minoris 


1.3 ±0.3' 


475 ± 30' 


4040 ± 100 K' 


0.7 ± 0.1 yuHz' 


18 + 9 days' 


Sun 


1 


1 


5780 K 


134 yuHz 


3.2 + 0.2 days' 



References: " Frandsen et al. (2002); * Stello et al. (2006); '' Bedding et al. (2006); '' CaiTier et al. (2007); ' Miglio & Montalban (2005); ' Bedding 
et al.(2004); « Fletcher et al. (2006); '' Kjeldsen et al. (2005); ' North et al. (2007); ' Bedding et al. (2007); * Appourchaux et al. (2008); ' Leccia 
et al. (2007); "' Arentoft et al. (2008); " Teixeira et al. (2008); ° Pijpers (2003); '' de Ridder et al. (2006); Barban et al. (2007); Tarrant et al. 
(2007); ' Tarrant et al. (2008); ' Chaplin et al. (2005). 




Fig. 1. Average mode lifetimes (t) from the stellar equilibrium and pulsation computations (diamonds) and from observations of 12 
stars (crosses), plotted as a function of: the L/M of each star (left-hand panel); and the effective temperature, Teff, of each star 



law. However, Houdek & Gough did make different assumptions 
about the "non-locality" of the convection to those adopted in 
the models used here. We are currently working on new models 
for a selection of red giant stars, which will provide a consistent 
treatment of the convection parameters. 

We therefore propose that a power law of the form 

(r) (X Tj 

does provide an adequate scaling relation for main-sequence, 
sub-giant and giant stars. 

Some of the measured lifetimes, close to T^ff ~ 5000 K, are 
noticeably shorter than the lifetimes implied by the best-fitting 
power law. These measurements were made by the Kjeldsen & 
Bedding group, and we suggest that the results for these stars 
may have been affected by not including rotation in the analysis 
(cf. Section|2|. 

The result we obtain when fitting the lifetimes from the stel- 
lar equilibrium and pulsation computations to a power law in 
L/M disagrees somewhat with the result of Chaplin et al. (2007). 
Chaplin et al. obtained A oc (L/M)' ^ while our result here 
implies A oc (L/M)"^. Chaplin et al. used pulsation computa- 
tions of a sample of models matching 22 solar-like stars in the 
Mount Wilson Ca ii H and K survey (Baliunas et al. 1995). All 
models had L/M lower than ~ 1.3-times the solar value. We 



suggest, with reference to Fig.[T] that this is a too small range in 
L/M to extract a robust relation. 

4. Discussion 

The aim of this study was to find a simple scaling relation be- 
tween the mean lifetimes, (r), of the most prominent solar-like p 
modes in stars, and the fundamental stellar parameters. We used 
predicted lifetimes from stellar equilibrium and pulsation com- 
putations, and measured mode lifetimes from observations of 12 
main-sequence, sub-giant and red-giant stars. Our study suggests 
that the lifetimes follow to a good approximation a scaling rela- 
tion that depends only on effective temperature; 

<T> ^ t;^. (3) 

The fourth-power dependence on effective temperature deserves 
some comment. The total surface flux radiated by a star is of 
course by definition proportional to T*^, and we would expect 
that, to first order, this will also correspond approximately to the 
surface convective heat flux in stars with near-surface convec- 
tion zones. The total heat flux transported by the convection is 
clearly an important parameter in fixing the strength of damping 
of solar-like oscillations. It is perhaps worth noting (see below) 
that a similar, simple dependence on effective temperature is not 
seen in the amplitudes of the solar-like p modes. 
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We may combine our scaling relation for average lifetime 
with the scaling relation of Kjeldsen & Bedding (1995) for the 
maximum mode amplitude. A, to yield a scaling relation for the 
maximum mode height, H, in the frequency power spectrum. For 
narrow-band intensity observations, Kjeldsen & Bedding pro- 
posed: 

L 

A = (6L/Lh ^ oc 



MT- 



ff 



(4) 



Combining with our prediction for lifetime, we obtain 

H (X g-\ 

We predict that the maximum mode height - which fixes the 
S/N, and hence the visibility, of peaks when they are resolved 
in the frequency power spectrum - depends predominantly on 
the surface gravity of stars, when the observations are made in 
intensity. We use the qualification in the preceding sentence be- 
cause the above does assume that the Kjeldsen & Bedding scal- 
ing for {6L/L)^ is accurate. While this scaling is supported to a 
large extent by observations, there is evidence for some depar- 
tures from the scaling, most notably for hotter solar-like stars (as 
noted in Section [1] above). When observations are instead made 
in Doppler velocity - which, using Kjeldsen & Bedding's pro- 
posed scaling, implies A oc L/M - we find 



Hoc 



rp4 

eff 



(5) 



Since the surface gravity changes fairly slowly along the 
main sequence, our new scaling relation for H given by 
EquationlH which assumes narrow-band intensity observations, 
suggests that stars notably cooler than the Sun might have mode 
heights that are in fact comparable to those of solar-like stars 
that are hotter than the Sun. Previously, asteroseismic target se- 
lections were based only on the expected amplitudes, mainly be- 
cause we had no idea about how the mode lifetimes scaled with 
fundamental stellar parameters. This would lead one to favour 
hotter stars as targets, while our result here suggests cooler solar- 
like stars may actually be very good targets as well. When obser- 
vations are instead made in Doppler velocity (Equation |5]l, our 
results imply that H increases with increasing effective temper- 
ature. 

Our prediction that cooler stars show longer lifetimes brings 
another added benefit: it is easier to estimate stellar rotation rates 
using asteroseismology, because the ratio between the rotational 
splittings and the mode linewidths will generally be higher than 
in hotter stars (see also Chaplin et al. 2008). This is important, 
because one of the main goals of asteroseismology is to mea- 
sure surface and internal differential rotation (see also: Gizon & 
Solanki 2003, 2004; and Ballot et al. 2006). 

We therefore suggest that in future our scaling relations for 
mode height (Equations |4] and |5]), and not the previously used 
relations for mode amplitude, be used when selecting targets for 
asteroseismology. 
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